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Cu ' The magnetic anisotropy af a rectangular Cu-O-Cu bond is investigated in second 

order of the spin-orbit interaction. Such a bond is characteristic for cuprates having edge 
sharing Cu02 chains, and exists also in the CU3O4 plane or in ladder compounds. For 
a ferromagnetic coupling between the copper spins an easy axis is found perpendicular 
to the copper oxygen plaquettes in agreement with the experimental spin structure of 
Li2Cu02. In addition, a pseudo-dipolar interaction is derived. Its estimation in the case 
of the CU3O4 plane (which is present for instance in Ba2Cu304Cl2 or Sr2Cu304Cl2) gives 

$-i( . a value which is however two orders of magnitude smaller than the usual dipole-dipole 

lyj ' interaction. 
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PACS. 75.30.Gw - Magnetic anisotropy. 
CN '■ PACS. 74.72.Jt - Other cuprates. 

^ ' PACS. 71.70.Ej - Spin-orbit coupling, Zeeman, Stark, and strain splitting. 
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To explain the rich magnetic behaviour of cuprates it is important to analyse not only the dominant 
isotropic superexchange interaction but also the smaller anisotropic terms due to the spin-orbit coupling 
0. In the case of a 180° Cu-O-Cu bond, this explains for instance the "easy plane" magnetism observed 
in the standard Cu02 plane of orthorhombic B or tetragonal M cuprates. However, there is another 
important structural element in many cuprates, namely a rectangular Cu-O-Cu bond. It is present in 
such cuprates having a CU3O4 plane (as in Ba2Cu304Cl2 Q or in Sr2Cu304Cl2 [||), having CUO2 chains 
(Li2Cu02 P,Rl, CuGeOs M) or in ladder compounds 0]. We are going to calculate here the parameters 
of the anisotropic superexchange interaction for the ideal 90° bond. To incorporate lattice effects we 
investigate a CU2O6 cluster (Fig. 1) - the common structural fragment of the cuprates mentioned. We 
choose as a representative example the Ba2Cu3 04Cl2 compound for which all important underlying 
parameters are derived from the extended multi-orbital tight-binding (TB) analysis of its bandstructure. 
Some results of the TB analysis were presented in a previous publication [Q . 

In the CU3O4 plane of A2CU3O4CI2 (A=Ba,Sr) one has to distinguish Cua with a local orthorhombic 
symmetry and Cub with local tetragonal symmetry. Lower than tetragonal local symmetry is obvious 
at the copper sites in the chains of Li2Cu02. This produces an additional feature - a weak higher order 
pseudodipole part of the magnetic anisotropy in these systems. In the leading orders of the perturbation 
theory the derived form of the anisotropic superexchange is similar for A2CU3O4CI2 and Li2Cu02. For 
CuGe03 the analysis should be completed by taking into account the influence of Ge side-groups g . 

The magnetic anisotropy can be calculated in perturbation theory. In second order of the spin-orbit 
coupling \(LaSa + LbSb) with a characteristic parameter A = O.leV for cuprates we obtain (see [pl|): 
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where J(myi, Ob; ^lAjOs) is the integral of superexchange interaction between the groundstate of Cub 
and the excited states tua and ua of Cua. The non-zero matrix elements of the orbital momentum are 
igj. = Lg = z, Lf^^ = —2i with the corresponding Hermitean conjugate ones. We will denote the 3d 
orbitals of Cu as follows: |0) = \x^ ~ y"^), |1) — |3z^ — r^), \z) — \xy), \x) = \yz) and \y) = \xz) (in the 
x-y coordinate system of Fig. 1). The small difference between the Cu^ and Cub crystal field splitting, 
which is inherent to the CU3O4 plane, has no meaning for the present purposes and we assume with 
high accuracy An ^ E^ - Eq ~ E^ - E^ and E^ ~ E^ = Eq. The values of A„ derived from the 
bandstructure data |10) of Ba2Cu304Cl2 are listed in Table 1 and these estimates are also relevant for 
other cuprates. 

To analyse (|l|,g) one has to calculate the different exchange integrals treating the Cu-0 hopping (~ tpd) 
and the 0-0 hopping (~ ipp,i7i-7r,ip7r,ipj) as perturbations. This is justified by the small ratios of these 
transfer integrals to the characteristic charge transfer energies (i?p,7r,p. — Eq). Together with the crystal 
field splitting the on-site interactions (involving direct and exchange Coulomb terms) at copper and 
oxygen ions are incorporated into the zero-order Hamiltonian. To avaluate the copper on-site Coulomb 
interactions the Racah parameters are taken as in Table 1. For B and C parameters the unscreened atomic 
values are used as explained elsewhere I14I . The parameter A is adopted in such a way that the Coulomb 
correlation C/oq_== A + 4i?-f 3C for the groundstate Cu-orbital coincides with the standard value (10.5 eV) 
given in Ref. ^^. Then one obtains also the other direct (C/^i = A - 45 -I- C, U^x,y = A-2B + C) and 
exchange [J^^ — AB + C, Jg^ = 3B + C) Coulomb integrals. For the correlation parameter Up at oxygen 
we take the standard value from Ref. |13] and Jp from m, which defines also Up^ = Up — 2Jp . The values 
of the transfer integrals are derived from the tight-binding analysis of the bandstructure flQ] . Finally, for 
proper estimates of the superexchange integrals one needs also the "bare" values for the charge transfer 
energies i?p,7r,p- ^ Eq, instead of screened ones inferred from the bandstructure calculations. To solve 



this problem we have used a cluster procedure |]T3 connecting self-consistently the mean-field (screened) 
values with "bare" ones for the on-site orbital energies. 

To perform calculations in a more effective way we introduce instead of the original in-plane oxygen 
\pi.2)~ ^-nd |7ri_2)-orbitals their combinations, like 
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with energies Ap/p = Ep—Eo^tpp and A^r/^ = ET^—EoTt-KTv, respectively (Table 1). For brevity we present 
below only the most important hybridization interaction giving rise to the lowest order contributions (4th 
order in tpd) to the superexchange 
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Contributions to the superexchange due to the oxygen orbitals at the right and the left cluster edges are 
mediated by 0-0 hybridization, which are taken into account below as the 5-th and 6-th order corrections. 

After these preliminaries we are able to estimate the different integrals J(m, 0;n, 0) entering into (^. 
Let us start, however, with the discussion of the basic isotropic superexchange between Cvla and Cus in 
the groundstate. In our notation that is given by J{0a,0b',0a,0b) = Jab- This exchange of 4th order 
in tpd and with a dominant ferromagnetic contribution due to the Hund's rule coupling Jp at oxygen was 
estimated already several times (see |g-|l^,|lj| ) . In the next order (6th) of perturbation theory we find 
an antiferromagnetic contribution which arises due to the hopping ip^ mn . With the parameters chosen 
that sums up to a resulting ferromagnetic exchange of Jab— -6 meV. 

Let us estimate now the other non-zero integrals responsible for the anisotropy ([3J3). That has to be 
considered as additional terms on top of the dominant isotropic ferromagnetic term discussed above. First, 
we calculate J{za,Ob', za,Ob) = J\ 
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Jazz with two contributions in accordance with two possible 



intermediate virtual states in the 4th order processes. These processes with the corresponding amplitudes 
and interaction parameters are shown schematically in Fig. 2. For the first contribution ~ Ja zz^ ^^^ 
singlet-triplet splitting in the virtual two-hole state, with occupied oxygen Jtt)- and |7r)-orbitals, is due to 
the correlation interaction Up. Then taking into account the other oxygen orbitals at the cluster edges, 
which gives the 5th order correction ~ tpp^ we obtain the result 
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with the triplet and singlet energies A^- = A^ -I- A#, Af- = A^ -I- A# -I- Up, respectively. The 5th order 
correction gives only a 10 per cent decrease. In the intermediate state of the second contribution J^ 1,^, 
the two holes meet at the same Cus ion, but in different orbitals \0b) and |1_b). They interact due to 
the Hund's coupling ~ Jq^ to give 
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without any correction of fifth order. The integral J{zb, Oa; zb, Qa) can be obtained by replacing tt, tt 
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p,p. Analogously we find J{xa,0b',xa,0b) = J a 
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where Ap^ = Ep^ — Eq is the on-site energy of the out-of-plane |p2)-orbital and £^ = A^ + 1^^. The same 
result also holds for J{xb,Oa] xb,Oa) and for J{yA, Ob; yA, Ob) = J{yB,OA] J/s, O^i) since Uq^ = U^y and 
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The pseudodipolar part of the superexchange for the Cu-O-Cu bond in a lattice arises due to deviations 
from the local tetragonal symmetry. Taking as a representative example the CU3O4 plane this deviation 
occurs at the Cu^ sites. The orthorhombic energy splitting between the \x')- and |j/')-orbitals (see the 
x'-y' coordinate system of Fig. 1) is derived from the bandstructure results at the F point to be Ssxy ~ 30 
meV. That corresponds to a mixed matrix element {x\H\y) = Ssxy/'i in the x-y coordinate system chosen 
in the calculation. The second source for the pseudodipolar interaction is the difference in the transfer 
integrals Stp_ = tp^ — tp^, between |p2)-orbitals perpendicular or parallel to the A-B bond line. A rough 
estimate derived from the bandstructure results |lG| is found to be Stp, « 40 meV. Using these deviations 
from tetragonality as an additional perturbation we obtain the following estimate for the only non-zero 
off-diagonal exchange integrals 

J{xA,OB;yA,OB) = J(yA,OB;xA,OB) ^ ^J^^\xa,0b;xa,0b) (9) 

starting from the 4th order term of the diagonal exchange ((0) and (H) without 5th order correction). 
The reduction factor is 

?=^ + 2|l^. (10) 
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Using the parameters for Ba2Cu304Cl2 (see Table 1) we obtain finally from (p|)-([l(]|) the following numbers 
for the anisotropic superexchange 

F^^ == -0.65meV, F^"= = F^^ == -10/ieV, jF'^^'l - 0.2/xeV . (11) 

In terms of the spin components parallel and perpendicular to the bond, the Hamiltonian (1^) takes the 
following form 

^ab = J\\SaSb + J±Sj^Sb + J^S^Sg , (12) 

with J|| = -F^^ -I- 2r='y, J_i = -F^^ - 2F^2' and .h = F^^ - 2F^^. It should be noted that in spite of the 
nontetragonality of the particular bond under consideration we did not find in the leading orders of the 
perturbation procedure any process which could lead to the Dzyaloshinskii-Moryia |11[| spin coupling. 

Let us discuss now the consequences for several cuprates and start with the CU3O4 plane. First, one 
has to compare the derived value of the pseudodipolar part of the anisotropic superexchange J*^ — 
(J|| — J_l)/2 — 2Txy with the estimate [5] Jdd ~ 20 /xeV for the usual dipole-dipole interaction for the 
Cu^-Cus bond. The small value of J^^ supports the suggestion |^ that the weak ferromagnetic moment 
observable in A2CU3O4CI2 is mainly due to the usual dipole-dipole interaction. Effects of the other kind 
of anisotropy given by AJ — (J|| -I- Ji^)l2 — J^ = — 2(F2z — T^x) > could be revealed only in a more 
general context. Actually, to predict the low-temperature magnetic structure in the above compound 
one has to take into account the main interations (including the isotropic and anisotropic terms) for 
Cu^-Cuyi and Cvlb-Cub bonds in the CU3O4 plane. This requires a more extended analysis for these 
two interacting antiferromagnetic subsystems ||,^. 

At the same time, in Li2Cu02, the estimate AJ > is important to explain the magnetic ordering. 
Based on preliminary bandstructure results |lq| one may infer that the relevant parameters for Li2Cu02 
will be only slightly changed in comparison to Tabic 1. Therefore, it can be expected that the above result 
for AJ remains qualitatively correct. At low temperatures, for a chain with ferromagnetically aligned 
spins, the anisotropy A J > prefers a spin direction parallel to the z axis, i.e. perpendicular to the chain 
plaquettes. That agrees with the experimental situation MM- The measured energy gap in the spin wave 



spectrum m leads to a uniaxial anisotropy which is in agreement with our estimate for AJ. But the 
values of the isotropic superexchange coupling reported inR] are about one order of magnitude smaller 
than that estimated from the bandstructure calculation ]16| and characteristic also for other cuprates. 
It would be interesting to reanalyse the spin wave spectrum starting from a spin Hamiltonian based on 
the microscopic electronic structure. Some additional complications should be noted as well. Actually, a 
magnetic frustration occurs in Li2Cu02 since the isotropic antifcrromagnetic superexchange for second 
neighbours within a chain is rather strong. That is expected |1™,16| to exceed even the ferromagnetic 
coupling for nearest neighbour spins. The stability of the overall magnetic structure in Li2Cu02 with 
antialigned ferromagnetic chains can be explained due to a considerable three dimensional spin coupling 
in this compound p^ , p6[ |. 
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TABLES 

TABLE 1: Parameters (in eV) to calculate the magnetic anisotropy. 
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FIGURES 



FIG. 1: The CU2O6 cluster used to calculate the magnetic anisotropy (full circles: copper, open circles: 
oxygen) with the ground state orbitals |0) at copper and the in-plane oxygen orbitals (full lines: \p)- 
orbitals and broken lines |7r)-orbitals). Also shown are the two coordinate systems mentioned in the 
text. 



FIG. 2: Examples for 4th order processes which contribute to jj, ' and J 
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